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Abstract. The Cloud Aerosol LIdar with Orthogonal Po- 
larization (CALIOP), on board the CALIPSO platform, has 
measured profiles of total attenuated backscatter coefficient 
(level 1 products) since June 2006. CALIOP’ s level 2 prod- 
ucts, such as the aerosol backscatter and extinction coeffi- 
cient profiles, are retrieved using a complex succession of 
automated algorithms. The goal of this study is to help iden- 
tify potential shortcomings in the CALIOP version 2 level 2 
aerosol extinction product and to illustrate some of the moti- 
vation for the changes that have been introduced in the next 
version of CALIOP data (version 3, released in June 2010). 
To help illustrate the potential factors contributing to the un- 
certainty of the CALIOP aerosol extinction retrieval, we fo- 
cus on a one-day, multi-instrument, multiplatform compari- 
son study during the CALIPSO and Twilight Zone (CATZ) 
validation campaign on 4 August 2007. On that day, we 
observe a consistency in the Aerosol Optical Depth (AOD) 
values recorded by four different instruments (i.e. space- 
borne MODerate Imaging Spectroradiometer, MODIS: 0.67 
and POLarization and Directionality of Earth’s Reflectances, 
POLDER: 0.58, airborne High Spectral Resolution Lidar, 
HSRL: 0.52 and ground-based AErosol RObotic NETwork, 
AERONET: 0.48 to 0.73) while CALIOP AOD is a factor 
of two lower (0.32 at 532 nm). This case study illustrates 
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the following potential sources of uncertainty in the CALIOP 
AOD: (i) CALIOP’s low signal-to-noise ratio (SNR) leading 
to the misclassification and/or lack of aerosol layer identifi- 
cation, especially close to the Earth’s surface; (ii) the cloud 
contamination of CALIOP version 2 aerosol backscatter and 
extinction profiles; (iii) potentially erroneous assumptions 
of the aerosol extinction-to-backscatter ratio (5 a ) used in 
CALIOP’s extinction retrievals; and (iv) calibration coeffi- 
cient biases in the CALIOP daytime attenuated backscatter 
coefficient profiles. The use of version 3 CALIOP extinction 
retrieval for our case study seems to partially fix factor (i) 
although the aerosol retrieved by CALIOP is still somewhat 
lower than the profile measured by HSRL; the cloud con- 
tamination (ii) appears to be corrected; no particular change 
is apparent in the observation-based CALIOP 5 a value (iii). 
Our case study also showed very little difference in version 
2 and version 3 CALIOP attenuated backscatter coefficient 
profiles, illustrating a minor change in the calibration scheme 
(iv). 


1 Introduction 

The Cloud Aerosol Lidar with Orthogonal Polarization 

(CALIOP), on board the CALIPSO platform (flying as part 

of the A-Train satellite constellation since April 2006), is a 

three-channel elastic backscatter lidar optimized for aerosol 
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and cloud profiling. CALIOP measures high-resolution 
(1/3 km in the horizontal and 30 m in the vertical at 532 nm 
in low and middle troposphere) profiles of the attenuated 
backscatter by aerosols and clouds at visible (532 nm) and 
near-infrared (1064 nm) wavelengths along with polarized 
backscatter in the visible channel (Winker et al., 2009). 
These data are distributed as part of the level 1 CALIOP 
products. The level 2 CALIOP products are derived from 
the level 1 measurements using a complex and intricate suc- 
cession of algorithms that are described in detail in a special 
issue of the Journal of Atmospheric and Oceanic Technology 
(e.g., Winker et al., 2009). The level 2 retrieval scheme is 
composed of a feature detection scheme, a module that clas- 
sifies features according to layer type (e.g., cloud vs. aerosol) 
and sub-type, and, finally, an extinction retrieval algorithm 
that estimates the aerosol backscatter, the extinction coeffi- 
cient profile and total column aerosol optical depth (AOD) 
for an “assumed” (see Sect. 2.2) extinction-to-backscatter ra- 
tio (also called S a ) for each detected aerosol layer. 

For a select list of observables, i.e., CALIOP attenuated 
backscatter, aerosol backscatter and extinction coefficient 
profiles have been shown to yield reasonable agreement with 
ground-based (Kim et al., 2008; Mamouri et al., 2009; Mona 
et al., 2009; Pappalardo et al., 2010) and airborne lidar mea- 
surements (McGill et al., 2007; Omar et al., 2009; Rogers et 
al., 2011). For example, Pappalardo et al. (2010) have ob- 
served a mean difference of less than 20 % between level 1 
CALIOP and ground-based EARLINET (European Aerosol 
Research Lidar Network) lidar measurements of attenuated 
backscatter profiles since June 2006 over Europe, showing 
an absence of evident biases in the CALIOP raw signals. 
Rogers et al. (2011) have conducted the most extensive quan- 
titative assessment study of the CALIOP 532 nm total attenu- 
ated backscatter to date, using coincident data from 86 under- 
flights by the NAS A-Langley High Spectral Resolution Lidar 
(HSRL) (Hair et al., 2008) acquired since June 2006. Results 
show HSRL and CALIOP (version 2) 532 nm total attenuated 
backscatter agree within 1.1% ±23% for daytime lighting 
conditions in the free troposphere. Kim et al. (2008) showed 
that CALIOP, when compared to a ground-based lidar in Ko- 
rea, has detected cloud and aerosol top/base layers and re- 
trieved the aerosol extinction profiles correctly within respec- 
tively 0.10 km and 30% in cloud- free nighttime and semi- 
transparent cirrus cloud conditions. According to Omar et 
al. (2009), CALIOP (version 2) generally overestimates the 
HSRL extinction measurements for several case studies, with 
an average extinction bias of 0.003 km -1 (~24%) during the 
CALIPSO and Twilight Zone (CATZ) validation campaign 
and 0.015 km -1 (~59%) during the Gulf of Mexico Atmo- 
spheric Composition and Climate Study (GoMACCS). 

Nonetheless, there are significant uncertainties associated 
with the version 2 CALIOP aerosol extinction and backscat- 
ter retrievals, and these are not well-quantified in any ancil- 
lary quality assurance information included in the level 2 data 
files. These uncertainties are introduced by several differ- 


ent factors that are often related to each other (Winker et al., 
2009; Yu et al., 2010). First of all, the CALIOP layer detec- 
tion scheme will most likely fail to detect layers with aerosol 
backscatter coefficients falling below a sensitivity threshold 
of 2~4xl0 4 km 1 sr 1 in the troposphere (Winker et al., 

2009) . Consequently, if we assume a lidar extinction-to- 
backscatter ratio (Sa) of 50 sr, the minimum detectable ex- 
tinction coefficient is on the order of 0.01 to 0.02 km -1 (cor- 
responding to a lowest detectable AOD of 0.02-0.04 in a ho- 
mogenous 2 km planetary boundary layer). A second sig- 
nificant source of error is the lack of photons returned from 
underneath highly attenuating layers, such as dense aerosol 
and cloud layers. This may result in the erroneous or total 
lack of aerosol identification in the lower part of a given pro- 
file. In such situations, the CALIOP detection algorithm can 
incorrectly identify the lower portions of an aerosol layer as 
being clear air, and thus no aerosol extinction coefficients 
will be reported for these regions. A third factor impact- 
ing the CALIOP extinction retrieval is the occasional mis- 
classification of layer type. Aerosols can be misclassified as 
clouds, and vice versa (Liu et al., 2009). Classification errors 
can also occur in the aerosol subtyping algorithm (Omar et 
al., 2009), leading to an incorrect assumption about the ap- 
propriate extinction-to-backscatter ratio. The CALIOP AOD 
fractional error is similar to the 5 a fractional error for small 
AOD values (Winker et al., 2009). On the other hand, as 
the AOD increases, the AOD fractional error will quickly be- 
come much higher than the 5 a fractional error. For example, 
a fractional error of 30% for S a would result in an AOD frac- 
tional error of ~50% for an AOD of 0.5 and nearly 100% for 
an AOD of 1 . 

Despite these uncertainties, there have been a number of 
publications using CALIOP version 2 level 2 data in a quali- 
tative or even quantitative manner. Focusing on articles pub- 
lished in 2010, some authors recognize the largely unvali- 
dated nature of level 2 version 2 data. Among those, there 
have been attempts to produce more accurate CALIOP data 
by applying further cloud-screening (Sekiyama et al., 2010) 
or even an intensive data screening scheme (Yu et al., 2010). 
Some mention the uncertainties associated with the level 2 
version 2 data but apply no specific filtering (e.g., Peyridieu 
et al., 2010, Jones and Christopher, 2010). We note that many 
articles in 2010 (and probably a few more in the previous 
years) omit discussions on the accuracy of level 2 version 2 
CALIOP data. This is, for example, the case for Gonzi and 
Palmer (2010), who qualitatively compared biomass burn- 
ing injection height estimates from the GEOS-Chem model 
to unfiltered CALIOP vertical feature mask data. This lat- 
ter product is also used to suggest the presence of an ex- 
tended aerosol layer over central India that could be associ- 
ated with agriculture crop residue burning activities (Sharma 
et al., 2010), and to help determine the altitude of smoke 
plumes over the US during Summer 2006 (McMillan et al., 

2010) . Finally, Kuhlmann and Quaas (2010) make more in- 
tensive use of the unfiltered level 2 CALIOP aerosol layer 
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product to draw conclusions regarding the particle type and 
general aerosol vertical distribution during the Asian sum- 
mer Monsoon. The conclusions drawn in these works are not 
necessarily wrong, and, in fact, may be absolutely correct. 
However, it is difficult to ascertain the merit of the reported 
results, simply because the data from which they are derived 
are not yet well validated. 

In this study we assess the consistency between the 
CALIOP AOD retrievals and comparative aerosol observa- 
tions from multiple sources and platforms (including ground- 
based, airborne and satellite instruments). We focus on a one- 
day, multi-sensor case study (Sect. 3) that was part of the nine 
ground-based CATZ field campaigns (each campaign occur- 
ring on separate days between 26 June and 29 August 2007) 
in Virginia and Maryland, when four AERONET sites were 
deployed and the NASA Langley Research Center airborne 
HSRL was flown along the daytime CALIOP track, with 
coincident space-borne observations available from MODIS 
and POLDER (Polarization and Directionality of Earth’s Re- 
flectances). The detailed suborbital observations, and in par- 
ticular, the comparison of coincident CALIOP and HSRL 
profiles, are used to explore the following potential rea- 
sons for the overall bias between the MODIS AOD and the 
CALIOP version 2 AOD product: (i) CALIOP’s low signal- 
to-noise ratio (SNR) which can lead to the misclassification 
and/or lack of aerosol layer identification, especially close to 
the Earth’s surface; (ii) the cloud contamination of CALIOP 
aerosol backscatter and extinction profiles; (iii) a potentially 
erroneous S& assumption in CALIOP’s extinction retrieval 
and (iv) calibration errors in the CALIOP daytime attenuated 
backscatter coefficient profiles. 

We note that this study is not a global validation assess- 
ment of the CALIOP version 2/version 3 aerosol extinction 
product but focuses on one case study. It is not intended 
to provide a quantitative assessment regarding the generality 
and magnitude of errors likely to be present in either version 
of the CALIOP extinction retrieval algorithm. Instead, our 
study identifies shortcomings in the version 2 level 2 aerosol 
extinction product and illustrates some of the motivation for 
the changes that have been introduced in the next version of 
CALIOP data (Version 3, released in June 2010). 

Based on the multi-instrument, multi -platform comparison 
study, we seek to quantify the major factors potentially con- 
tributing to the uncertainty of the CALIOP aerosol extinction 
retrieval. We then use version 3 CALIOP measurements to 
assess a potential correction of each of these factors. We sub- 
mit that the identification and discussion of retrieval uncer- 
tainties provided here will help understand and interpret the 
results obtained in previous studies like the ones cited above. 


2 Instruments 

2.1 AERONET 

The AErosol RObotic NETwork (AERONET) (Holben et 
al., 1998) is composed of automatic sun- sky scanning spec- 
tral radiometers that determine AOD by direct sun measure- 
ments. AERONET measurements also provide the Angstrom 
exponent (A), which expresses the wavelength dependence, 
X, of the AOD and is defined as the slope of the first order lin- 
ear regression of log(AOD) versus log(k). The aerosol size 
distribution and optical parameters (such as the single scat- 
tering albedo, volume concentration, refractive index, etc.) 
are derived from the angular distribution of sky radiances 
measured in the almucantar according to the algorithm de- 
veloped by Dubovik and King (2000a). In this study, we use 
version 2-level 1.5 AERONET data (Smirnov et al., 2002). 
During the CATZ field experiment, the AERONET sunpho- 
tometer observations were sampled more frequently than in 
the case of the standard automatic mode measurement pro- 
tocol (Holben et al., 1998), preventing the data from being 
labeled level 2. However, the correct calibration of the sun- 
photometers during the experiment results in the same esti- 
mated total uncertainty in the direct AOD measurements as 
for the level 2 data: ~0.010-0.021 (Eck et al., 1999). In 
the case of AOD values above 0.2 at 440 nm, Dubovik et 
al. (2000b) reports accuracies of 0.03 for the single scatter- 
ing albedo, 0.02-0.04 for the real part of the refractive index, 
30% (50%) of the imaginary part of the refractive index in 
case of low (high) absorption, 15-35% (15-100%) of the vol- 
ume size distribution in case of a radius between 0. 1 and 7 pm 
(lower than 0.1 pm or above 7 pm). In the case of lower AOD 
values (AOD(440)<0.2), the accuracy levels drop down to 
0.05-0.07 for SSA, 80-100% for the imaginary part of the 
refractive index, and 0.05 for the real part of the refractive 
index. 

2.2 CALIOP 

CALIOP on the CALIPSO platform employs a linearly po- 
larized laser that transmits pulses at 532 nm and 1064 nm. 
The two 532 nm receiver channels separately measure the 
components of the 532 nm backscatter signal polarized par- 
allel and perpendicular to the outgoing beam. The measured 
CALIOP attenuated backscatter coefficient at wavelength X 
and range r, /^(r), can written as: 

p[(r) = (fi a , x (r) + p m ,x(r))T^(r) (1) 

where and are, respectively, the aerosol and molec- 
ular backscatter coefficient profile, and T x = T^ x 7j^ k T ^ x 
is the atmospheric two-way transmittance (i.e., signal atten- 
uation) due to aerosols, molecular scattering, and absorbing 
gases such as ozone. 
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Table 1 . CALIOP parameters used in this study with attributed name, variable, file, level and horizontal/vertical resolution- version 2 in black 
and version 3 in red italic. The parameters /? 532 @ 1/3 km, cloud @ 1/3 km, aerosol @5 km and feature @5 km have the same horizontal and 
vertical resolution in CALIOP version 2 and version 3. 


Name 

Variable 

File 


Level Resolution due to averaging 

Horizontal Vertical ( < 8 km) 

p' S32 @1/3 km 

Total _Attenuated_Backscatter_532 

CAL_LID_L1 -Prow (ValStagel )-V2(V3)-01.. .hdf 


1 

1/3 km 

30 m 

cloud @1/3 km 

Layer _Top .Altitude ; Layer _B ase .Altitude 

C AL_LID_L2 .3 3 3mCLay-Prov- V2( V3 )-0 1. . .hdf 


2 

1/3 km 

30 m 

aerosol @5 km 

Layer .Top Altitude ; Layer _B ase .Altitude 

CAL_LID_L2_05 kmALay-Prov-V2(V3)-01 ...hdf 


2 

5 km 

30 m 

feature @5 km 

Feature _Clas sification .Flag s 

CALXIDX2_VFM-Prov(Va/Stage7)-V2(V3)-0L..hdf 


2 

5 km 

30 m 

Pa , 532 @40 km (@5km) 

Total_B ackscatter _Coefficient_5 32 

C ALXID _L2 .40 km AProCal -B eta- V2( 5 kmA Pro-Prov- V3)-0\ . . 

.-.hdf 

2 

40 (5) km 

120 (60) m 

a a , 532 @40 km ( @5 km) 

Extinction .Coefficient .5 3 2 

CAL_LID_L2_40 kmAProCal-Beta-V2(5 kmAPro-Prov- V3)-01. . 

.-.hdf 

2 

40 (5) km 

120 (60) m 

S a> 532 @40 km (@5km) 

“a a ,532 @40f5jkm”/“^ a 532 @40(5) km” 

None 


2 

40(5) km 

120 (60) m 


The two-way aerosol transmittance between the range of 
the lidar calibration region, r c , and range r, T^ x (r c ,r), can 
be expressed as follows: 

T lx 6 c,r) = ex P ^-2 J oc a ,x(r )dr j =exp(- 2 r a ,0 ( 2 ) 

where a a ,A is the aerosol extinction coefficient and r a ,i, the 
aerosol optical depth. 

We will be concentrating mostly on the CALIOP- 
measured total attenuated backscatter coefficients, /3 532 (r), 
the retrieved aerosol backscatter coefficients, /l a , 532(4), and 
the retrieved aerosol extinction coefficient profiles, a a? 532(4) 
along the CALIOP ground track at 532 nm. 

The extinction coefficient profiles are retrieved using a 
globally automated feature recognition algorithm that as- 
sumes a range-invariant extinction-to-backscatter ratio, also 
referred to as lidar ratio (^532 = a a? 532 O')/ At, 532(7*)) for each 
layer detected. The CALIOP value of 6^532 used for any 
layer depends on the geographical location, the integrated at- 
tenuated backscatter color ratio, the layer-integrated volume 
depolarization ratio, and a general look up table (LUT) (Liu 
et al., 2009; Omar et al., 2009). The “assumed” CALIOP 
lidar ratio in this study is, in fact, selected after several 
steps including the subtype classification of each detected 
aerosol layer, based itself on an extensive clustering analy- 
sis of global AERONET measurements. The prelaunch goal 
of the CALIPSO mission was to retrieve aerosol extinction 
coefficients accurate to within ±40% (Winker et al., 2003). 
We have attributed names to all the CALIOP parameters used 
in this study. They are listed in Table 1 along with the stan- 
dard variables, original file name, level, and spatial resolution 
due to averaging for CALIOP ’s version 2 and version 3 data 
products. We emphasize that this study is mainly focused 
on the CALIOP version 2 extinction products. CALIOP ver- 
sion 3 data products are used specifically to address whether 
the new CALIOP algorithm fixes certain issues identified in 
version 2. 

CALIOP’ s version 2 data products do not provide un- 
certainty estimates for retrieved optical parameters such as 
AOD and extinction coefficients (see CALIPSO Version 2 
Data Quality Statements). The uncertainties attributed to the 


CALIOP aerosol optical depths can be obtained by apply- 
ing an error estimator algorithm to the quantities reported in 
the aerosol layer products, taking into account the relative er- 
ror on the lidar ratio, the calibration coefficient and the SNR 
for each detected aerosol layer. The error in the extinction 
retrieval due to the SNR may be slightly more complex to 
estimate as it depends on the backscatter intensity, the light- 
ing conditions (i.e., day vs. night), and the amount of hori- 
zontal averaging applied to the initial attenuated backscatter 
profiles. 

2.3 HSRL 

Retrieval of aerosol extinction profiles using the standard 
elastic backscatter lidar technique requires either a measure- 
ment of AOD to constrain the extinction retrieval (Young, 
1995; McGill et al., 2003) or an assumption on the aerosol 
extinction-to-backscatter ratio value (Cattrall et al., 2005). 
On the other hand, the HSRL technique directly measures 
the vertical profile of aerosol extinction and extinction-to- 
backscatter ratio, without requiring ancillary aerosol mea- 
surements or assumptions about aerosol type (Hair et al., 
2008). The HSRL technique is typically employed for the 
532 nm wavelength utilizing the iodine vapor filter technique 
(Hair et al., 2001, 2008; Piironen and Eloranta 1994). The 
received 532 nm backscatter return is split between two op- 
tical channels to discriminate between aerosol and molec- 
ular returns: (1) one passing through an iodine vapor cell 
which absorbs the central portion of the backscatter spec- 
trum, including all of the Mie backscatter, and transmits 
only the Doppler/pressure-broadened molecular backscatter 
(called the “molecular channel”); and (2) another one that 
passes all frequencies of the signal returned by both aerosols 
and molecules (called the “total scatter channel”). In addi- 
tion, the lidar is polarization-sensitive, enabling discrimina- 
tion between spherical and non- spherical particles. The first 
channel (“molecular channel”) is used to retrieve the extinc- 
tion profile and both channels are used to retrieve profiles 
of aerosol backscatter and extinction coefficients and aerosol 
depolarization ratios. Hair et al. (2008) described the poten- 
tial errors introduced in any of these quantities and found the 
532 nm extinction systematic errors to be less than 0.01 km -1 
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Table 2. HSRL parameters used in this study with attributed name, variable and short description. 


HSRL data (subset file: ~4/3 km horizontal and 30 m vertical resolution) 

Name 

Variable 


Description (at 532 nm) 

/ 

^532 

532 .total _attn _b sc _cloud_screened 

Attenuated backscatter coefficient with cloud mask applied 
Reported from ~60 m to ~1.5 km below the plane (~7.5 km) 

£a,532 

532 _bsc _cloud_screened 

Aerosol volume backscatter coefficient with cloud mask applied 
Reported from ~60 m to ~0.5 km below the plane (~8.5 km) 

«a,532 

532_ext 


Retrieved 532nm aerosol extinction coefficient. 

Reported from ~360 m to ~2.5 km below the plane (~6.4 km) 

S a,532 

Sa_532 


Retrieved Extinction-to-backscatter ratio 

Reported from ~360 m to ~2.5 km below the plane (~6.4 km) 

AOD 532 

AOTJii 


Aerosol optical thickness determined from the molecular 

channel near the aircraft and near the surface 

Derived from ~60 m to ~2.5 km below the aircraft (~6.4 km) 


for typical aerosol loading. Table 2 describes the HSRL an- 
alyzed data products used in this study. We use an HSRL 
subset file with a ~4/3 km horizontal and 30 m vertical reso- 
lution. On the one hand, the ~4/3 km horizontal resolution of 
the HSRL aerosol backscatter coefficient profiles is obtained 
by computing 10 s running averages of the raw data (initially 
sampled at 2 Hz), then sub-sampling the results by a factor 
of 20. On the other, the ~4/3km horizontal resolution of 
the HSRL extinction and lidar ratio coefficient profiles is ob- 
tained by computing 60 second running averages of the raw 
data. 

2.4 POLDER and MODIS 

POLDER-3 (POLarization and Directionality of Earth’s Re- 
flectances, the third version of the POLDER instruments, on 
board the PARASOL platform) and MODIS (on board the 
Earth Observing System (EOS) AQUA satellite) are both 
passive radiometers, with both platforms being part of the 
A-Train during five years (December 2004-2009), includ- 
ing our study period of Summer 2007. POLDER’s strength 
is the measurement of directional and polarized characteris- 
tics of the solar radiation reflected by the Earth- Atmosphere 
system. MODIS, on the other hand, has a finer spatial 
and spectral resolution. POLDER AOD estimates of po- 
larizing particles over land surfaces are retrieved in the 
865 nm channel. MODIS AOD is retrieved over oceans in 
7 different spectral bands (6+1 extrapolated) from the vis- 
ible to the near infrared and over land in 3 bands (2+1 
interpolated). POLDER’s spatial resolution is 5x6.5 km 
(500 x500 m for most MODIS channels) and its wide field 
of view induces a 1600 km swath (2330 km for MODIS) that 
allows a nearly global daily coverage. To increase the sig- 
nal to noise ratio, the standard retrieval algorithm is applied 
to 3x3 POLDER pixels (20x20 for MODIS), leading to a 


resolution in the aerosol AOD of 15x19.5 km (10x10 km 
for MODIS at nadir). The AOD retrieval from the POLDER 
polarized measurements is described by Deuze et al. (2001) 
and the MODIS AOD retrieval algorithm over land is de- 
scribed in Kaufmann and Tanre (1998). The polarization by 
aerosols mainly comes from small spherical particles in the 
accumulation mode (Vermeulen et al., 2000), indicating that 
POLDER-derived AOD is well suited for remote sensing of 
fine mode particles. Validation studies suggest that the ex- 
pected uncertainty on the MODIS AOD over dark land sur- 
faces could be represented by AAOD = ±0.05ib0.15 AOD 
(Levy et al., 2010). 


3 Evaluation of Version 2 CALIOP extinction retrieval: 
4 August 2007 (A CATZ case study) 

3.1 Aerosol type and sources 

The MODIS true color RGB image in Fig. la shows some 
haze hovering over a significant part of the Mid- Atlantic East 
Coast of the United States, extending from Virginia to New 
Jersey on 4 August 2007. This particle plume is most likely a 
mix of aerosol pollution from regional anthropogenic sources 
and smoke coming from wildfires in the Northwestern United 
States. According to the National Interagency Fire Center, 
more than a dozen large fires were reported from late July 
to early August of 2007 in the Northern Rockies of Idaho 
and Montana. By 07 August, those fires had affected nearly 
400 000 acres in Idaho and had produced smoke that blan- 
keted much of the United States. The 3 day-HYSPLIT air 
mass back-trajectories at three different heights from 500 to 
1500m (Draxler et al., 2010; Rolph, 2010) (Fig. la), sug- 
gests that a part of the aerosol plume over the East Coast on 
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Table 3. 4 August 2007 - Distance (km) between the different measurements during the CATZ experiment; from top to bottom, the CATZ- 
Sanders AERONET station (white diamond, Fig. 2a), the closest point on the HSRL track to CATZ-Sanders (green line, Fig. 2a), the closest 
point on the CALIOP ground track to CATZ-Sanders (white line, Fig. 2a), the closest CMAQ cell to CATZ-Sanders containing a MODIS 
and CALIOP version 2 AOD value (red box, Fig. 2a) and the closest CMAQ cell to CATZ-Sanders containing a POLDER AOD retrieval 
(black box, Fig. 2a). 


Distance (km) 

CATZ.Sanders Closest point on 
HSRL track 

Closest point on 
CALIOP track 

CMAQ cell with 
CALIOP/MODIS AOD 

CATZ .Sanders 

- 

- 

- 

- 

Closest point on HSRL track 

0.940 

- 

- 

- 

Closest point on CALIOP track 

0.138 

0.908 

- 

- 

CMAQ cell with CALIOP/MODIS AOD 

5.809 

5.315 

5.680 

- 

CMAQ cell with POLDER AOD 

17.703 

17.339 

17.569 

12.067 


4 August 2007 may have come from the Northern part of the 
United States. 

We will be focusing our study over the CATZ-Sanders El- 
ementary School AERONET station (39.04° N; -77.51° W), 
one of the four sunphotometer sites that were deployed along 
the CALIOP track during the CATZ campaign. This station, 
shown by a white diamond in Fig. la, will hereafter be called 
“CATZ-Sanders”. CATZ-Sanders was 138 m away from the 
CALIOP track and the overpass on 4 August occurred at 
18:27 UTC. Aerosol microphysical and optical properties de- 
rived from the inversion of two angular sky-radiance mea- 
surements at CATZ-Sanders on 4 August 2007 are shown in 
Fig. lb. 

The aerosol plume over CATZ-Sanders seems predom- 
inantly composed of fine particles, with Angstrom coeffi- 
cients (A between 440-870 nm) of 1.92 (Fig. lb). This is 
confirmed by the volume size distributions that show, for 
both measurements, a peak around 0.16 pm in radius. Fi- 
nally, the particles show significant light absorption with a 
single scattering albedo coefficient (<z>o) between 0.94 and 
0.96 and an imaginary part of the refractive index (Im(^)) of 
about 0.01. 

3.2 Ground-based, airborne and space-borne AOD 
measurements 

In this study, for convenience, all satellite data are remapped 
on the 12x12 km Community Multiscale Air Quality 
(CMAQ) model grid (US EPA, 1999). Each MODIS 
10 xlOkm cell center has been attributed to the closest 
CMAQ cell center. In the case of CALIOP, the product to 
be remapped is the standard level 2 extinction coefficient, 
a a ,532 @40 km (see Table 1). CALIOP provides one con- 
stant extinction vertical profile between start-location l start 
and end-location l en d> with a horizontal distance of 40 km be- 
tween l start and lend- A 12 x 12 km CMAQ cell can then con- 
tain, at the most, two different parts of 0^532 @40 km pro- 
files. When the CMAQ cell contains only one a df 532 @40 km 
profile, this profile is simply attributed to the cell’s cen- 
ter; on the other hand, when the cell contains two differ- 


ent a a> 532 @40 km profiles, the final profile attributed to the 
cell’s center is an average of those two 0^532 @40 km pro- 
files weighted by the number of corresponding /3 532 @ 1/3 km 
profiles contained in the cell. The CALIOP AOD data value 
for each cell is then obtained by integrating its corresponding 
a a? 532 @40 km profile vertically. 

Figure 2a shows the locations of CATZ-Sanders (white di- 
amond), the CALIOP ground track along the closest 40 km 
segment (white line), the corresponding airborne HSRL 
track segment (green line) and the closest CMAQ 12 x 12 km 
cell (red box). Recall that all satellite data are remapped 
onto the CMAQ grid and the closest CMAQ cell to CATZ- 
Sanders (red box in Fig. 2a) contains a remapped MODIS 
and CALIOP AOD observation. On the other hand, the clos- 
est CMAQ cell with available POLDER AOD data on 4 Au- 
gust 2007 is shown as a black box in Fig. 2a, at a distance of 
~18km between CATZ-Sanders and the closest POLDER 
extinction observation. Table 3 summarizes the horizontal 
distances between each measurement during the experiment. 

Figure 2b shows the collocated ground-based (sunpho- 
tometer, black), airborne (HSRL, orange) and space-borne 
(MODIS green, POLDER red and CALIOP blue) AOD ob- 
servations. The sunphotometer is the only instrument pro- 
viding a full temporal evolution of AOD values throughout 
the afternoon and evening of 4 August 2007. The HSRL in- 
strument completes this temporal information with two over- 
passes over CATZ-Sanders around 16:48 and 17:52 UTC. 

It should be noted that HSRL overflew CATZ-Sanders 
~30min before CALIPSO overpass time (17:52 UTC com- 
pared to 18:27 UTC for CALIPSO) and ~900m away from 
the CALIOP ground- track (Table 3). A ground-based wind 
profiler instrument in Belts ville (Maryland) shows an aver- 
age wind speed of ~2.6m per second from the surface up 
to ~3.8 km between the HSRL and CALIOP overpass time, 
blowing mainly from the northwest. Accordingly, a 30-min 
air mass travel time between the HSRL and CALIOP ob- 
servations would represent a distance of roughly 5 km at the 
ground. 
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Fig. 1. 4 August 2007 - (a) MODIS true color RGB image over the East Coast of the United States, 72 h-HYSPLIT air mass back- trajectories 
at 500 (red), 1000 (blue) and 1500 (green) meters above model ground level arriving at 19:00 UTC over CATZ-Sanders (white diamond) and 
the CALIOP ground track (white line); the back- trajectories are computed using the gridded meteorological data archives of the National 
Weather Service’s (NWS) National Centers for Environmental Prediction (NCEP) Global Data Assimilation System (GDAS) model; (b) 
Version 2-level 1.5 CATZ-Sanders AERONET observations derived from the angular distribution of sky radiance in the almucantar (2 
measurements). 


Whether it is statistically relevant to compare AOD re- 
trievals from HSRL and CALIOP on that time (~30 nm) and 
horizontal scale (~5 km) is difficult to ascertain. The auto- 
correlation of either the HSRL AOD or the CALIOP AOD 
dataset for different lagged distances along the flight or the 
satellite track could help inform on a relevant comparison 
distance. For CALIOP and HSRL, the distances of interest 
to this study (a few tens of kilometers) would likely be sam- 
pled in a short enough time that it would be reasonable to 
neglect temporal evolution of the aerosols. On another hand, 
the statistical relevance of comparing HSRL and CALIOP on 
a 30-min time scale could be partially informed by the au- 


tocorrelation coefficient of the CATZ-Sanders ground-based 
AERONET AOD values for different lagged times. 

Instead of performing the studies described above, we re- 
fer to (Anderson et al., 2003) that uses datasets from four 
different instruments (two ground-based, one airborne and 
one space-borne) to show a reasonably coherent picture of 
mesoscale aerosol variations: when measured at one local 
point, scales below 4.4 h or 70 km show large autocorrelation 
coefficients ( R above 0.90), indicating both small ambient 
variability and measurement noise. Redemann et al. (2006) 
shows similar results ( R =0.96) on a scale of 15 km for 
the AOD derived from the NASA Ames Airborne Tracking 
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Sunphotometer (AATS-14) during the Extended-MODIS-A 
Validation Experiment (EVE) campaign in April of 2004. 

The relevance of a 30 min and 5 km spatio-temporal scale 
between CALIOP and HSRL AOD also varies from one en- 
vironment to another. In our case study, according to Fig. 2b, 
there is a fair amount of variation in the AERONET AOD 
measurements throughout the afternoon and evening of 4 
August 2007 (from 0.48 to 0.87 at 532 nm). The variation 
zb 1/2 h around the time of the A-Train overpass is smaller 
but still significant, with AOD values (at 532 nm) ranging 
from 0.48 to 0.73. This variation, similar to the range of 
AOD 1/2 h preceding the A-Train overpass, corresponds to a 
change of ~35% in the AOD (0.25 compared to 0.71 at the 
A-train overpass time) over a course of ~5 km (distance cov- 
ered by the air mass in 1/2 h with an averaged wind speed 
of ~2.6m/s). The temporal variability of midvisible AOD 
during our case study of 4 August 2007 falls in the AATS- 
14 spatial AOD variability range (2-37%) observed along 
most 6 km- segments through an environment dominated by 
biomass-burning plumes over Canada during the ARCTAS 
field campaign (Arctic Research of the Composition of the 
Troposphere from Aircraft and Satellites) (Shinozuka et al., 
2010). In addition, the effect of vertical mixing could either 
decrease or increase variability in remotely sensed total col- 
umn AOD observations as well as profiles of aerosol extinc- 
tion properties. 

The HSRL AOD retrieval (0.52) agrees with the 
AERONET direct-beam AOD measurement (0.57), showing 
a slight difference of 0.05 at the time of the second HSRL 
overpass (~18:00UTC). 

At the time of the A-train overpass (dashed vertical line 
in Fig. 2b), MODIS and AERONET report similar AOD re- 
trieval values (with a difference of 0.04, below MODIS’s 
AOD uncertainty of ~0.15, 0.05 + 15% of 0.67). On the 
other hand, POLDER underestimates the AERONET AOD 
by 0.13. This slight difference could be due to uncertainties 
in the POLDER inversion algorithm. Some bias could also be 
due to the satellite’s coarse spatial resolution in a temporally 
and spatially varying aerosol field, especially for POLDER 
with a coarser resolution than MODIS and further away from 
the sunphotometer (Table 3). Recall that POLDER is sensi- 
tive to fine polarizing particles over land and, thus, retrieves 
the fine mode AOD when MODIS retrieves the AOD corre- 
sponding to the entire volume size distribution of the parti- 
cles (see Fig. lb). 

In conclusion, all three AOD observations [i.e., MODIS 
(0.67), POLDER (0.58) and HSRL (0.52)] are contained in 
the AERONET AOD envelope within 1/2 h around the A- 
Train overpass (0.48 to 0.73 at 532 nm). The CALIOP ver- 
sion 2 AOD value (0.32), however, is well below the range 
of all other AOD measurements in Fig. 2b. 

A broader comparison between version 2-derived 
CALIOP AOD and collocated MODIS AOD retrievals (both 
on the CMAQ grid) over the Continental United States 
during the summer of 2007 yielded (results not shown 


here) a similar mean underestimation by CALIOP. The 
best statistical agreement between those two quantities 
was found over the Eastern part of the United States with, 
nonetheless, a weak correlation (R~0.4) and an apparent 
CALIOP version 2 underestimation (by ~66 %) of MODIS 
AOD. 

In the following section, we investigate the potential rea- 
sons for a disagreement between the AOD calculated from 
CALIOP ’s version 2 extinction product and the rest of the 
AOD measurements in Fig. 2b. 

3.3 HSRL and CALIOP backscatter and extinction 
coefficient profiles 

Figure 3 shows the CALIOP and HSRL P 522 cross sections 
of attenuated backscatter (also called “curtain scene”) along 
the 40 km segment of their ground tracks close to CATZ- 
Sanders (respectively corresponding to the white and green 
lines in Fig. 2a). Both CALIOP and HSRL are shown at a 
horizontal resolution of ~4/3km (output resolution of the 
“subset” HSRL file and sliding average of four CALIOP 
P 532 @1/3 km profiles). A dashed black vertical line in all 
three panels shows the closest profile to CATZ- Sanders on 4 
August 2007. 

The difference between the CALIOP “curtain scenes” 
shown in Fig. 3a and b reflects an additional cloud- screening 
of the data. Yost et al. (2008) compared MODIS images 
overlaid with the CALIOP cloud @ 1/3 km product, and the 
feature @5 km product. It was shown that the CALIOP 1/3- 
km detection results are entirely consistent with the MODIS 
image. However, in regions populated by broken bound- 
ary layer clouds, layers detected at coarser resolutions (1- 
km and above) are frequently misclassified as cloud. This 
was determined to be strictly a coding error in the cloud- 
clearing procedure, and not related to the algorithm design. 
To circumvent this error, in this part of our study, an ad- 
ditional cloud screening has been applied to all CALIOP 
ft 532 @ 1/3 km profiles using the cloud @ 1/3 km product: all 
CALIOP P 532 @1/3 km coefficients are deleted underneath 
the highest detected cloud in the cloud @1/3 km product. 

The black circle in Fig. 3a, b and c points out a region of 
the curtain scene showing high initial raw fi 532 @ 1/3 km co- 
efficient values around 2.2 km on the vertical (Fig. 3a). This 
signal is classified as a cloud in the cloud @ 1/3 km prod- 
uct and is removed in Fig. 3b, thanks to the additional cloud 
screening described above. In Fig. 3b, the sliding average of 
four profiles in the curtain scene (for an end horizontal reso- 
lution of 4/3 km) prevents us from observing a lack of data 
underneath the detected cloud, where all P' 522 @ 1/3 km coef- 
ficients were removed. Figure 3c shows a lack of HSRL data 
in the corresponding region, most probably due to the pres- 
ence of clouds as well (the HSRL data are cloud- screened, 
see Table 2). 

Figure 3 illustrates the differences in the SNR of the 
HSRL and CALIOP instruments. The CALIOP curtain scene 
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Fig. 2. 4 August 2007 - (a) Co-localization of CATZ-Sanders (white diamond), the closest 40km-CALIOP track segment (white line), 
the corresponding closest 40 km-HSRL track segment (green line), the 12 x 12 km closest CMAQ cell (red box) reporting available MODIS 
and CALIOP version 2 AOD data and the closest CMAQ cell reporting available POLDER AOD data (black box); (b) Temporal evolution 
of AERONET CATZ-Sanders direct sun AOD measurements (black) and the collocated HSRL AOD retrievals during the two overpasses 
(AOD 532 of Table 2, orange); At the time of the A-Train overpass (18:27 UTC, dotted grey line) over CATZ-Sanders: in green, the MODIS 
AOD retrieval (red box of (a)), in red, the POLDER AOD value (black box of (a)) and in blue, the CALIOP version 2 AOD value (red box 
of (a)). All AOD observations are either retrieved or computed at 532 nm (use of the Angstrom exponent between 440-675 nm, 470-670 nm 
and 865-670 nm for AERONET, MODIS and POLDER). 


(Fig. 3b) appears much noisier than the HSRL cross section 
(Fig. 3c), which makes it harder to analyze in terms of po- 
tential atmospheric vertical composition. On the other hand, 
Fig. 3c seems to show two fairly separate and spatially homo- 
geneous stronger regions in the /3 532 intensity on the vertical: 
the lowest one lies roughly between 1 and 2 km and the up- 
permost one is around 3 km. In addition, the closest point 
to CATZ-Sanders on the HSRL track (black dashed line in 
Fig. 3c) seems fairly representative of the rest of the 40 km 
curtain scene. 


Figure 4a2 shows the most closely collocated version 2 
CALIOP and HSRL fi' 532 profile to CATZ-Sanders (black 
dashed line in Fig. 3b and c). Both CALIOP (Fig. 4a2, 
blue) and HSRL (Fig. 4a2, red) profiles are shown at a 
~4/3km resolution (output resolution of the HSRL sub- 
set file and selection of the closest CALIOP profile in the 
4/3 km-resolution curtain scene of Fig. 3b). The CALIOP 
/3 532 profile still clearly shows a very low SNR compared to 
the HSRL p' 532 profile. 
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Fig. 3. 4 August 2007 - CALIOP version 2 raw (a), CALIOP ver- 
sion 2 with further cloud- screening (b) and HSRL (c) ^ 532 curta i n 
scene along the 40 km segment of its ground tracks close to CATZ- 
Sanders (respectively white and green line in Fig. 2a). Both 
CALIOP and HSRL are shown at a ~4/3 km resolution; the black 
circles in (a), (b) and (c) point out a region with cloud contamina- 
tion. 


CALIOP’s low SNR (as shown in Figs. 3b and 4a), es- 
pecially in daytime, requires spatial averaging of the attenu- 
ated backscatter profile over a significant horizontal distance 
to detect potential features. This is one of the tasks of the 
Selective Iterated Boundary Locator (SIBYL) in CALIOP’s 
automated level 2 product routine (Vaughan et al., 2009). In 
short, SIBYL consists of an algorithm that iteratively aver- 
ages profiles at different horizontal scales (5, 20 or 80 km), 
scans those averaged profiles to detect aerosol and cloud lay- 
ers, and removes detected layers from the profiles before fur- 
ther averaging. As a result, strongly scattering layers and 
portions of layers are detected at finer spatial resolution, 
while more tenuous regions are detected at coarser resolu- 
tions. All layers detected are then classified according to type 
and subtype (Liu et al., 2009; Omar et al., 2009). Particulate 
backscatter and extinction coefficients are then derived for 
each layer detected at the 5-km, 20-km, and 80-km averag- 
ing interval, using profiles of (z) averaged horizontally to 
the spatial resolution at which the layer was detected (Young 
and Vaughan, 2009). In CALIPSO’s version 2 data products, 
the level 2 “native resolution” yd a , 532 and a a ,532 profiles are 
further averaged (layers detected at 5-km or 20-km) or repli- 
cated (80-km layers) as required to be reported at a uniform 
final resolution of 40 km horizontal and 120 m vertical (Ta- 
ble 1). 

The closest version 2 ^,532 @40 km profile to CATZ- 
Sanders is shown in Fig. 4b2 (black), along with the col- 
located HSRL /3 a ,532 profile (red). Unlike the processing 
of CALIOP profiles, we saw no necessity to average the 
HSRL profiles on a similar horizontal distance at the ground 
because of HSRL’s considerably higher SNR and accuracy. 
Figure 3c supports this decision by showing a spatially uni- 
form atmospheric HSRL curtain scene in the vicinity of 
CATZ-Sanders. In addition, the HSRL would cover 40 km 
in a few minutes (HSRL flies at ~1 17 m/s) compared to a 
few seconds for CALIOP (flies at ~7 km/s), adding potential 
temporal differences in the HSRL-CALIOP comparison. 

In Fig. 4b2, the CALIOP /3 a? 532 @40 km profile reports 
no aerosol above ~3.2km or below ~ 1.4 km. Both the 
CALIOP /J a> 532 @40 km and the HSRL /3 a? 532 profile seem 
to show mostly two intensity peaks in the vertical. The 
change in intensity between the uppermost and the lowest 
peak could be due to either a change in the particle type (size 
and shape, hence different aerosol cross section and phase 
function) and/or a change in the particle concentration, and 
does not necessarily show two vertically separated aerosol 
layers. Concerning the uppermost aerosol peak, the HSRL 
and CALIOP signals compare fairly well between 2.3 and 
3.2 km. In Table 4, the standard aerosol @5 km product (Ta- 
ble 1), shows that the uppermost aerosol layer in version 2 
(detected at a horizontal averaging of 20 km) is located be- 
tween 2.7 and 3.1 km, and defines it as polluted dust aerosol 
particles (CALIOP model 5^ = 65 sr). The lowest intensity 
peak in Fig. 4b2 consists of a fairly constant portion of the 
HSRL /3 a? 532 profile recording roughly 0.003 km -1 sr -1 from 
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Fig. 4. 4 August 2007 - (a2) HSRL (red) and CALIOP version 2 (blue) /^ 32 profile, both at ~4/3 km resolution, (b2) HSRL ^532 (red) and 
CALIOP version 2 ^ 532 @40 km profile (black), c2) HSRL 5 a? 532 (red) and CALIOP version 2 5^532 @40 km profile (black) and (d2) 
HSRL o ' a? 532 (red) and CALIOP version 2 0^532 @40 km profile (black). All profiles are the closest to CATZ-Sanders. The CALIOP and the 
HSRL track are ~900 m away from each other and the HSRL overpass is ~30 min early (17:52 UTC) compared to CALIOP (18:27 UTC); 
(a3) through (d3) are the same as (a2) through (d2) but using version 3 for all CALIOP profiles. 
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the lowest few hundred meters close to the ground up to 
1.9 km. Although the corresponding CALIOP profile starts 
around 1.4 km and misses a lot of the aerosol signal ob- 
served by the HSRL, it seems to pick up the lowest peak with 
an overestimation of 1 x 10 -3 km -1 sr -1 at 1.9 km before 
a maximum of 5.9xl0 -3 sr -1 km -1 at 2.2km. According 
to Table 4, the standard CALIOP version 2 aerosol @5 km 
products define the lowest aerosol layer (detected with a hor- 
izontal averaging of 80 km) as being located between 1.5 and 
2.5 km, and composed of dust aerosol particles (CALIOP 
model S d = 40 sr). 

The presence of either polluted dust or pure dust aerosol 
particles as detected in the version 2 CALIOP retrieval algo- 
rithm is unlikely according to the findings of Sect. 3.1. In- 
deed, the optical and microphysical properties of the aerosol 
plume over CATZ- Sanders tend to show a predominance of 
fine and strong light absorbing particles, possibly coming 
from a mix of haze and biomass burning particles. 

In summary, Fig. 4b2 shows fairly good agree- 
ment between the HSRL /3 a , 532 and CALIOP version 2 
/3 a ,532 @40 km profiles, except for a lack of CALIOP val- 
ues below ~ 1.4 km and a strong peak in the CALIOP ^,532 
signal around 2.2 km. The immediate reasons could be that 
(i) CALIOP, with its low SNR, cannot detect tenuous aerosol 
layers or reach all the way down to the lidar-detected sur- 
face due to aerosol attenuation and (ii) there is a significant 
bug in the cloud-screening algorithm, that could explain the 
disparity between CALIOP and HSRL ^,532 around 2.2 km 
(corresponding to the height at which a cloud is reported in 
Fig. 3a). 

Figure 4c2 compares the CALIOP version 2 6^532 
@40 km profile (= 0^532 @40km//3 a? 532 @40 km, see Ta- 
ble 1) with the measured HSRL £^532 profile (see Table 2). 
For HSRL, 5^ 532 b) is simply the ratio of a d , 532 (z) and 
^a,532 (z), where both quantities are measured directly by the 
instrument. The random uncertainty for the HSRL S a values 
are determined using the same methodology as presented in 
Hair et al. (2008). For the profile in plot 4c2, the uncertainty 
in S a is less than 20 % (above 1 km) except at the very top 
of the profile (>3 km) where the aerosol loading approaches 
zero and hence relative errors are large. The increased un- 
certainties in this case study are somewhat larger than those 
presented previously (Hair et al., 2008) due to the relatively 
high AOD (0.52) and the presence of clouds in the sampling 
interval. Applying the requisite cloud clearing for this case 
leads to fewer raw HSRL profiles being averaged, which in 
turn reduces the SNR of the averaged HSRL S a profile (see 
Sect. 2.3). As a consequence, the HSRL S d variability in 
Fig. 4c2 should not be attributed solely to the natural vari- 
ability of aerosols. The CALIOP retrieval algorithm does not 
assume a profile of S a values but assumes, instead, a single 
S d value for each detected aerosol layer in the vertical. The 
fact that the 6^532 @40 km profile in Fig. 4c2 varies in the 
vertical is due to the averaging of aerosols that were iden- 
tified as different types and detected at different horizontal 


scales. Although CALIOP and HSRL show similar averaged 
S d values in the vertical (66 sr for CALIOP compared with 
64 sr for the HSRL), CALIOP shows a much smaller range of 
S a> 532 @40 km (from 56 to 70 sr) compared with the HSRL 
(from 53 to 83 sr between 3.2 and 1.4 km). The reason is that 
the range of different S d values attainable in the CALIOP au- 
tomated algorithm is much smaller than in reality. This ob- 
servation leads to the introduction of a third potential expla- 
nation for the discrepancies between CALIOP and the HSRL 
extinction observations: (iii) the “assumed” CALIOP ^,532 
value for each aerosol layer detected in the vertical could be 
erroneous and shows less variability than in reality. 

The small variation of the CALIOP 6^532 @40 km profile 
in Fig. 4c2 explains the strong resemblance of the CALIOP 
Ai, 532 @40 km and 0^532 @40 km profiles in Fig. 4b2 and 
d2. The HSRL 0^532 profile in Fig. 4d2 clearly shows an in- 
crease in the extinction coefficient values between 2.4 and 
3 km, followed by a stronger peak extending from ~2 km 
down to a few hundred meters close to the ground. On the 
other hand, the CALIOP a ^532 @40 km profile reports the 
uppermost increase higher than for the HSRL with an ap- 
proximate difference of 500 m in the vertical and seems to 
pick up ~500m of the lowest aerosol peak (between 1.4 and 
1.9 km). 

To summarize, there are several important dissimilarities 
between the CALIOP version 2 and the HSRL extinction co- 
efficient profiles on 4 August 2007. The potential reasons for 
those discrepancies are investigated in the remainder of this 
study. We also investigate whether these factors have been 
addressed in the version 3 CALIOP extinction product. 

3.3.1 CALIOP ’s detection of tenuous aerosol layers and 
signal attenuation 

We attempt to estimate the impact of failed detection of 
low-level aerosol layers due to high signal attenuation on 
column AOD, using the collocated HSRL a ^332 profile of 
Fig. 4d2 (red). The integration of the HSRL 0^532 profile 
from the ground to the base of the lowest layer detected 
by CALIOP (leading to an AOD of 0.23 from a few hun- 
dred meters to 1.5 km), and again beginning above the top 
of the highest layer detected by CALIOP (AOD of 0.01 
from 3 km to the top) adds a total of 0.24 to the standard 
CALIOP AOD of 0.32. This would, at least, account for the 
amount of extinction needed for CALIOP to be consistent 
with the AERONET AOD range 1/2 h around the overpass 
(0.48 to 0.73) on 4 August 2007. Based on the strong MODIS 
AOD underestimation (~66 %) by version 2-CALIOP de- 
rived AOD we have observed (results not shown here) over 
the continental Eastern United States during the summer of 
2007, we speculate that this is not just a problem specific to 
our case study but instead, may occur on larger space and 
time scales. Indeed, the CALIOP team has developed an al- 
ternative retrieval philosophy for low-lying aerosol layers. In 
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Table 4. 4 August 2007 - Version 2 and version 3 CALIOP 5 km-aerosol layer products (layer top and base altitude, averaging required to 
detect each layer, layer feature sub-type classification and quality assurance) for the closest 5 km-profile to CATZ-Sanders along the CALIOP 
track. 


CALIOP aerosol @5 km 

version 2 

version 3 

[Layer Top-Layer Base] (km) 

[3. 1-2.7] [2.5-1.5] 

[3. 1-2.7] [2.5-1. 5] [1. 3-0.2] 

Averaging for layer detection (km) 

80 20 80 

80 20 80 

Layer feature sub-type classification 

dust polluted-dust 

dust polluted-dust polluted-dust 

Layer feature sub-type QA 

high high 

high high low 

Initial layer lidar ratio (532 nm, sr) 

40 65 

40 55 55 


those cases where transparent aerosol layers are detected, if 
(a) the initial estimate of layer base is “close to” the Earth’s 
surface, and (b) the surface is reliably detected, and (c) the 
mean attenuated backscatter between the initial base estimate 
and the surface is positive, the layer base estimate is revised 
downward to a new, lower altitude very near the surface. This 
new scheme has been implemented in version 3 data prod- 
ucts, and preliminary results suggest that it will have the de- 
sired effects (Vaughan et al., 2010). 

Concerning the version 3 CALIOP algorithm, the very first 
result relevant to this section is the presence of aerosols be- 
low 1.9 km in Fig. 4b3, c3 and d3, compared to no detec- 
tion of the lowest aerosol layer for CALIOP version 2 in 
Fig. 4b2, c2 and d2. This is confirmed in Table 4 by an ad- 
ditional aerosol layer close to the Earth’s surface (between 
1.3 and 0.2km) in the CALIOP version 3 aerosol @5 km 
products. From 1.9 km down to —700 m, CALIOP version 
3 Pa, 532 @5 km (Fig. 4b3) stays within zblxlO -3 km -1 sr -1 
of the HSRL /? a ,532 profile, with the exception of a strong sig- 
nal of 5.4 xlO 3 km 1 sr 1 at — 1.5 km. Underneath 1.4 km, 
CALIOP version 3 Pa, 532 @5 km underestimates the HSRL 
Pa, 532 profile by up to — 2x 10 -3 km -1 sr -1 at an altitude of 
—200 m. We observe an overall better agreement between 
CALIOP version 3 ^532 @5 km and HSRL 0^532 (Fig. 4d3) 
compared with version 2 below 1.9 km (Fig. 4d2). 

However, the integration of <2^532 @5 km in the vertical 
yields a total column version 3 CALIOP AOD of 0.32, equal 
to the version 2 CALIOP AOD in Fig. 2b (0.32). In other 
words, applying the new layer base determination scheme 
used in the version 3 CALIOP retrievals does not induce any 
change in the total column AOD for our case study. The 
lack of improvement in the version 3 CALIOP AOD esti- 
mate (relative to version 2) is due to offsetting effects from 
other changes between version 2 and version 3, which are 
discussed in more detail below. 

3.3.2 CALIOP’s “assumed” lidar extinction-to- 
backscatter ratio value per detected aerosol 
layer 

An alternative CALIOP version 2 0^532 @40 km* profile 
was computed by applying a newly devised extinction re- 


trieval to all previously cloud-screened CALIOP version 2 
P532 @ 1^3 km profiles in the 40 km region of interest (such 
as shown in Fig. 3b with a —4/3 km horizontal resolution). 
The alternative extinction retrieval uses a simple iterative nu- 
merical method, starting from a range ro (here, correspond- 
ing to an altitude of —4 km) down to the ground. The aerosol 
extinction coefficient is assumed equal to zero at range ro, 
the molecular extinction and backscatter coefficient profiles 
are taken from the GEOS-5 model provided in the CALIOP 
level 1 data, and the 6^532 profile is taken from the closest 
HSRL profile to CATZ-Sanders (Fig. 4c2, red). Additional 
information on the alternative extinction retrieval is given 
in Appendix A. The alternative CALIOP AOD values along 
the 40 km segment are then obtained by integrating each al- 
ternative extinction coefficient profile in the curtain scene 
between —1.4km and —3.2 km, i.e., the range of CALIOP 
detected aerosol layers and extent of the standard CALIOP 
<*a ,532 @40 km profile in Fig. 4d2 (black). The result is a 
40 km-averaged alternative CALIOP version 2 AOD value 
of 0.44 compared to the standard CALIOP AOD value of 
0.32 close to CATZ-Sanders on 4 August 2007 (Fig. 2b). It 
appears that, in this case study, modifying the extinction-to- 
backscatter ratio profile in the CALIOP extinction retrieval 
has less of an effect on the final AOD retrieval (adds 0.12 
in the AOD) than the impact of failed detection of low-level 
aerosol layers due to high signal attenuation (adds 0.27 in the 
AOD, previous section). The conclusion of a minor impact 
of CALIOP’s potentially erroneous “assumed” S a value com- 
pared to the inability of its signal to reach all the way down to 
the surface on the AOD retrieval cannot yet be stated in gen- 
eral. This result may, indeed, be strongly influenced by very 
similar averaged HSRL and CALIOP S a values (Fig. 4c2) on 
4 August 2007 close to CATZ-Sanders. 

Table 4 shows no change in the version 3 aerosol sub-type 
classification of the two uppermost aerosol layers (dust and 
polluted dust). The layer close to the Earth’s surface is also 
classified as polluted dust. Among the changes in version 
3, the extinction solver uses a different “observation-based” 
lidar ratio S a for polluted dust: 55 sr at 532 nm instead of 
65 sr for version 2. Paradoxically, the use of this lower lidar 
ratio in the version 3 processing negates the improvements 
in layer optical depth that would have been expected from 
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Fig. 5. 4 August 2007 -HSRL / 3^ 2 (red), alternative CALIOP 
version 2 non-cloud- screened P^ncs @40 km* (magenta), al- 
ternative CALIOP version 2 cloud- screened p^ 2 cs @40 km* 
(green), and alternative CALIOP version 2 cloud- screened normal- 
ized ^532, cs @40 km** profile (black); p'£ 2 ncs @40 km* , ^ 2 cs 
@40 km* and P^ 2 cs @40 km** are obtained by a sliding average 
of four P 522 @ 1/3 km profiles followed by averaging all valid pro- 
files on the 40 km segment (white line in Fig. 2a); P^ 2 cs @40 km* 
and p^ 2 cs @40 km** are cloud- screened using the cloud @ 1/3 km 
product; p^ 2 cs @40 km** is normalized by the ratio of the mean 

HSRL p^ 2 and the mean p^ 2 cs @40 km* from 4.5 to 7.5 km in 
altitude. All profiles are collocated with CATZ-Sanders. 

the implementation of the revised base determination scheme 
(previous section). 

Version 3 CALIOP S a ,532 @5 km (Fig. 4c3) shows less 
variability (numerous repetition of two values in the vertical: 
40 and 55 sr) than version 2 (Fig. 4c2) due to less horizon- 
tal averaging (different horizontal resolution in Table 1). In 
addition, the vertical mean version 3 CALIOP 6^532 @5 km 
(~52 sr) agrees less well with the mean HSRL S a ,532 value 
(~64 sr) than the averaged version 2 ^,532 @40 km (~66 sr). 
In the end, version 2 and 3 seem to show a similar classifi- 
cation of the aerosol type and underestimation of the corre- 
sponding lidar ratio for our case study (Table 4). The latter 
contributes to the underestimation of the HSRL 0^,532 profile 
by the CALIOP version 3 a a ,532 @5 km profile in Fig. 4d3. 

3.3.3 CALIOP’s cloud clearing, averaging and 
calibration of the attenuated backscatter 
coefficient profile 

Figure 5 shows the closest HSRL p 522 profile (red) to 
CATZ-Sanders on 4 August 2007, along with three alterna- 


tive CALIOP version 2 /^ 32 profiles. The first one, called 
^532 ncs @40 km* (in magenta in Fig. 5), is obtained by ap- 
plying a sliding average of four P 522 @ 1/3 km profiles be- 
fore averaging all valid profiles in the 40 km segment close to 
CATZ-Sanders (white line in Fig. 2a). The second one, called 
^532 cs @40 km* (in green in Fig. 5), corresponds spatially to 
the first one, but with a sliding average of four profiles on the 
cloud-screened fi 522 @1/3 km curtain scene (Fig. 3b). 

We note that the first two alternative CALIOP profiles 
(magenta and green, Fig. 5) show more general variability 
than the HSRL profile (red, Fig. 5), illustrating the dif- 
ferences in SNR between the two instruments, and emphasiz- 
ing the utility of using a broader horizontal averaging scale of 
80 km as the input of CALIOP’s standard multi-scale averag- 
ing engine and feature detection algorithm. In addition, the 
comparison between CALIOP ^532 ncs @40 km* (magenta) 

and /^3 2 cs @40 km* (green) confirms the presence of a re- 
ported cloud in the 40 km of interest around 2.2 km in height. 

Concerning version 3, a closer look at Fig, 4b3 shows 
a strong decrease (~1.4x 10 -3 km -1 sr -1 ) in the inten- 
sity of CALIOP /3 a ,532 @5 km around 2.2km compared 
to the version 2 ^,532 @40 km peak (corresponding to 
5.9 x 10 -3 km -1 sr -1 in Fig. 4b) that we presumed to be due 
to cloud contamination. 

Two major factors need to be considered when compar- 
ing the HSRL p" 2 (red) and the CALIOP p' 5 c 22 cs @40 km* 
(green) profiles in Fig. 5. First of all, the instruments differ 
regarding their calibration technique and accuracy. The ac- 
curacy of the CALIOP level 1 products (and, by extension, 
many of the level 2 products) critically depends on the accu- 
racy of the calibration of the attenuated backscatter profiles. 
The nighttime CALIOP 532 nm parallel attenuated backscat- 
ter measurement is calibrated by determining the ratio be- 
tween the measured signal and the total backscatter estimated 
from an atmospheric scattering model (Powell et al., 2009; 
Hostetler et al., 2006; Russell et al., 1979) across an altitude 
range of 30-34 km, where aerosol loading is assumed to be 
low and there is still sufficient molecular backscatter to pro- 
duce a robust signal. Because of the degradation of the SNR 
in the calibration region due to noise associated with the so- 
lar background signal, the CALIOP daytime 532 nm calibra- 
tion coefficients are interpolated from the adjacent nighttime 
data segments (Powell et al., 2010). On the other hand, the 
Airborne HSRL is internally calibrated to a high accuracy 
(~ 1—2%), and does not rely on normalization to estimated 
backscatter from assumed clear-air regions for calibration 
(Hair et al., 2008). 

Secondly, the HSRL fi^ 2 (red) and CALIOP 

^532 cs @40 km* (green) profiles differ in terms of the 
atmospheric attenuation of each lidar signal. The attenuation 
of the CALIOP profile is measured relative to the base 
of CALIOP’s molecular normalization region at 30-km 
(the minimal beam attenuation above this region is in- 
cluded in the calibration coefficient). Because the HSRL 
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is internally calibrated, and does not rely on molecular 
normalization, atmospheric attenuation of the HSRL signal 
is measured relative to a point 1.5 km below the aircraft, 
zhsrl (~7.5 km). 

For those cases where there are no clouds above the HSRL, 
the magnitudes of the attenuated backscatter profiles mea- 
sured by the two instruments will differ by a factor of 

AT 2 = exp( —2 / (a m (z')+ao 3 (z'))dz ) (3) 

V ./30km / 

so that 

/J532, cs @40km*(z) = A T 2 pg 2 (z) (4) 

Aerosol loading is considered negligible between 30-km and 
Zhsrl, and thus no aerosol attenuation term is included in 
the calculation of AT 2 . The requisite values for <203 (z) and 
a m (z) are estimated using the gridded ozone and molecular 
number density profile data from the GEOS -5 analysis prod- 
uct available from the NASA Goddard Global Modeling and 
Assimilation Office (GMAO). 

AT 2 for the ^532 cs @40 km* profile in Fig. 5 is 0.88 
(molecular and ozone optical depth are respectively ~0.04 
and ~0.02). Hence, if the CALIOP signal was correctly cal- 
ibrated, HSRL (red) would be ~12% higher than the 
CALIOP /353 2 cs @40 km* (green) profile. Figure 5 shows, 

in fact, a general overestimation of the HSRL /3^ 2 profile 
(red), especially along the uppermost and lowest intensity 
peak. We have computed the difference between the inte- 
grated red and green profiles of Fig. 5 as follows: 


(/ z S H r^24 Z -/gr m ^2,cs @40km*<fe) x 100 


rZHSRL O'H d 
Jz=0km^532 az 


= 13.74% (5) 


The amount of overestimation of the integrated HSRL /3^ 2 
on the integrated CALIOP ^532 cs @40 km* profile is sim- 
ilar to what would be expected in the case of a correctly 
calibrated CALIOP signal on 4 August 2007 near CATZ- 
Sanders. Nonetheless, (Powell et al., 2010 and Rogers et al., 
2011) show that, in general, the CALIOP calibration remains 
an issue in the version 2 level 1 attenuated backscatter prod- 
ucts. Indeed, (Powell et al., 2010) showed that the use of a 
constant scaling factor to transfer calibration from nighttime 
to daytime measurements in the version 2 data products was 
precluded by thermally-induced misalignment of the trans- 
mitter and receiver, causing the daytime signal levels to vary 
non-linearly. In conclusion, the sub-optimal version 2 day- 
time calibration of CALIOP ’s raw signal can be added as a 
fourth potential reason for discrepancies between CALIOP 
and either HSRL or any other available AOD measurement. 
The next version 3 CALIOP data release improves upon this 
calibration scheme with significant modifications (Powell et 
al., 2010). 

For our case study, the version 3 CALIOP cloud-screened 
closest P 532 profile to CATZ-Sanders in Fig. 4a3 shows (in 


blue) very little difference to the corresponding version 2 
P 532 profile (Fig. 4a2) on 4 August 2007. The absolute mean 
difference between both version 2 and version 3 profiles is 
3.4x 10 -5 km -1 sr -1 (corresponding to ~4% of the mean 
version 3 /3 532 ). The absolute maximum difference between 
both blue profiles of Fig. 4a2 and a3 is 1.3 x 10 -4 km -1 sr -1 
(~14% of the mean version 3 P 532 )- 

As a result of the different signal attenuation, in order 
to rigorously compare the CALIOP and HSRL total atten- 
uated backscatter coefficients we normalized the CALIOP 
profile (yd 5 3 2 cs @40 km**, in black in Fig. 5) using the ra- 
tio of the mean /3^ 2 by the mean fi^ 2 cs @40 km* in a “clear 
air” region (from 4.5 to 7.5 km). The normalized CALIOP 
^532 cs @40 km** (black, Fig. 5) is fairly close to the HSRL 
P^ 2 profile (red, Fig. 5) with ~93% of the differences be- 
tween both profiles less than 0.5 x 10 -3 km -1 sr -1 . The inte- 
gration of both profiles on the vertical is within 1 % of each 
other. 

We note that the normalized CALIOP p'^ 2 cs @40 km** 

(black) and the HSRL /3^ 2 profile (red) should show the 
same trend, given that both instruments are sampling the 
same aerosol layer at the same wavelength. Both profiles 
decrease rapidly with altitude at heights below ~ 1.9 km, 
most probably due to strong aerosol attenuation. The 
HSRL backscatter and extinction profiles corresponding 
to the /3^ 2 profile (red) are measured directly while the 
CALIOP backscatter and extinction profiles corresponding 
to ^532 cs @40 km** (black) are only retrieved in those re- 
gions where an aerosol layer is identified. The strong aerosol 
attenuation of the signal below ~ 1.9 km in Fig. 5, together 
with the additional noise that CALIOP has to content with, is 
what causes the version 2 CALIOP layer detection algorithm 
to fail to identify the full vertical extent of the layer. This 
leads to a premature CALIOP apparent aerosol base height 
determination explaining the lack of aerosol reported below 
~ 1.4 km in Fig. 4b2. 

Figure 6 shows the comparison of the collocated HSRL 
(red, Fig. 5) and the renormalized CALIOP (black, Fig. 5) 
total attenuated backscatter coefficients from ~8 km down to 
the surface close to CATZ-Sanders on 4 August 2007. 

According to Fig. 6, a large number of collocated HSRL 
f3^ 2 and CALIOP fi^ 2 cs @40 km** coefficients are below 
1.5xl0 -3 km -1 sr -1 (57%). CALIOP shows a fairly simi- 
lar amount of those lower values of HSRL total attenuated 
backscatter coefficients below and above the one-to-one line 
with a comparable mean value of 8-9x10 4 km 1 sr l . On 
the other hand, the overall regression line in Fig. 6 (red 
line, first principal component regression method (Kendall, 
1957)) shows a slight CALIOP overestimation of the HSRL 
total attenuated backscatter coefficients. Nonetheless, after 
constraining the averaging, the cloud screening, and the nor- 
malization, the CALIOP level 1 attenuated backscatter mea- 
surements show good agreement with the HSRL /3^ 2 profile 
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(correlation coefficient of 0.91, insignificant offset and a 
slope very close to 1 ). 

4 Conclusion 

While first attempting to assess the general consistency be- 
tween both space-borne CALIOP column integrated aerosol 
extinction profiles and MODIS AOD retrievals, we observed 
low correlation (R~0.4) and a general underestimation (by 
66 %) of the MODIS-derived AOD by CALIOP (version 2) 
during Summer 2007 over the Eastern part of the United 
States. Both data sets are subject to retrieval error. The 
possible reasons for the discrepancies between both satellite 
retrievals as they pertain to CALIOP are discussed and ex- 
plored based on a carefully selected case study containing 
detailed multi-sensor, multi-platform aerosol observations 
(ground-based AERONET sunphotometer, airborne HSRL 
lidar and spaceborne POLDER/PARASOL, MODIS/AQUA 
and CALIOP/CALIPSO). The case study is on 4 August 
2007 and part of the CATZ field campaign over Maryland. 
On that day, we observed a consistency in the AOD values 
recorded by MODIS (0.67), POLDER (0.58), HSRL (0.52) 
and AERONET (0.48 to 0.73) while CALIOP was a factor 
of two lower (0.32 at 532 nm), hence our focus on the po- 
tential shortcomings of the CALIOP data product. The study 
of the coincident HSRL and CALIOP profiles helps illustrate 
what are likely to be the most important reasons for the over- 
all bias in the CALIOP version 2 AOD. Here we summarize 
these reasons and state whether they have been addressed in 
CALIOP version 3: 

(i) CALIOP’s low SNR (as shown in Figs. 3b and 4a) pre- 
vents the detection of tenuous aerosol layers. Further- 
more, as shown above, the attenuation of the signal by 
dense aerosol plumes can drive the signal within a layer 
below CALIOP’s detection threshold, and thus prevents 
identification of the full vertical extent of the layer. This 
explains the lack of CALIOP version 2 /3 a ,532 @40 km 
data below ~ 1.4 km and the premature determination of 
the aerosol layer base in our case study. Using the collo- 
cated HSRL layer AOD above and below the CALIOP 
detected aerosol layer altitudes adds a total of 0.27 to 
the CALIOP AOD value. CALIOP version 3 uses an 
alternative retrieval for low-lying aerosol layers. It re- 
vises the layer base estimate downward for certain con- 
ditions. In our case study, version 3 adds an aerosol 
layer close to the Earth’s surface, from 1.3 km down to 
200m. Both CALIOP version 3 backscatter and extinc- 
tion coefficient profiles show a better agreement with 
the corresponding HSRL profiles. However, because 
the CALIOP lidar ratio is still underestimated by ~20% 
(55 sr vs. ~70 sr for HSRL), the CALIOP version 3 re- 
trievals still underestimate the HSRL aerosol extinction 
and backscatter coefficients below 1 .4 km. 



Fig. 6. 4 August 2007 - HSRL / 3^ 2 ( re d profile in Fig. 5) 

versus alternative cloud- screened normalized CALIOP version 2 
^532 cs ® km** (black profile in Fig. 5) coefficients from ~8 km 
down to the surface close to CATZ-Sanders; First principal com- 
ponent regression method (red line): 532 , cs = (L 19=1=0.03) 

^532 + (0.00±0.00), R = 0.91, RMSD = 0.34 x 10~ 3 , N = 240. 

(ii) The natural variability of S a ,532 is not properly 
represented by the “observation-based” CALIOP 
S a> 532 @40 km profile used in the version 2 aerosol re- 
trieval algorithm. Applying an alternative extinction re- 
trieval to the version 2 CALIOP attenuated backscatter 
profiles using the independently measured HSRL lidar 
ratio profile adds ~0.1 to the total CALIOP AOD. The 
CALIOP version 3 6^,532 @5 km profile still shows very 
little variation and underestimates even more the coin- 
cident HSRL S a profile for our case study. The latter 
results from a misclassification of the aerosol type as 
desert dust (Initial CALIOP model S a of 40 sr for both 
version 2 and 3 at 532 nm) and polluted-dust (CALIOP 
model S a of 65 sr for version 2 and 55 sr for version 
3 at 532 nm) instead of what seems to be a predomi- 
nance of fine and strongly light absorbing particles, pos- 
sibly coming from a mix of haze and biomass burning 
aerosols (with a CALIOP model S a value of 70 sr for 
version 2 and 3). We would like to stress the potential 
importance of this factor. Further investigation and vali- 
dation of the CALIOP “assumed” 6^532 product should 
be carried out on a broader scale and time period (i.e. 
measurements of 6^532 along the CALIOP track over a 
large seasonal and spatial range). 

(iii) There is a significant bug in the version 2 cloud- 
screening algorithm. This presumably explains the 
disparity between CALIOP and HSRL measured total 
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attenuated backscatter coefficient, and in consequence, 
retrieved aerosol backscatter and extinction coefficient 
profiles around 2.2 km near CATZ-Sanders on 4 August 
2007. The cloud contamination has the opposite effect 
of artificially increasing the AOD value in the general 
underestimation of the CALIOP AOD. This cloud con- 
tamination bug of CALIOP level 2 parameters is fixed 
in version 3 for our case study. 

(iv) Finally, although the CALIOP signal seemed to be fairly 
well calibrated during our case study, it is important to 
mention that the version 2 CALIOP daytime calibra- 
tion scheme has proven to be suboptimal, leading to 
bias errors in the 532 nm total attenuated backscatter 
and, in consequence, will propagate into the CALIOP 
aerosol extinction products (Rogers et al., 2011; Powell 
et al., 2010). Our case study shows very little differ- 
ence between version 2 and version 3 CALIOP attenu- 
ated backscatter coefficient profiles. 

Let us mention that multiple scattering, which is assumed 
to be negligible in the CALIOP level 2 aerosol algorithms 
could also be a potential reason for the extinction retrieval 
uncertainties. Indeed, multiple scattering effects are more 
significant in the case of spaceborne than airborne lidar sys- 
tems due to a larger footprint. They can alter the apparent 
extinction or transmittance of the medium, lead to depolar- 
ization of the returned signal, and can produce stretching of 
the return pulse. Nevertheless, the effects of multiple scat- 
tering seem to apply mostly in the case of dense dust plumes 
recording high AOD values. Based on Winker et al. (2003), 
in the case of aerosols other than large dust particles, multiple 
scattering is likely to contribute, at best a 10% uncertainty 
to the retrieval of aerosol extinction profiles. The error in- 
troduced by ignoring multiple scattering effects is negligible 
compared to a fractional error of 30% in the lidar S a ratio (re- 
sulting in an AOD fractional error of ~50% when the AOD 
is around 0.5). In the case of fresh, dense dust layers close 
to the source region, the analysis of airborne in situ size dis- 
tribution observations during SAMUM-1 (Saharan Mineral 
Dust Experiment, Southern Morocco, May-June 2006) have 
shown that the multiple-scattering-related underestimation of 
the extinction coefficient in the CALIOP lidar signals ranges 
from 10-40% (Wandinger et al., 2010). On the other hand, 
Liu et al. (2010) shows that for moderately dense dust cases 
(AOD^l and extinction smaller than 1km -1 ), the vertical 
homogeneity of the particulate depolarization ratio profile in- 
dicates negligible impact from multiple scattering. Multiple 
scattering effects are not considered in our paper as all the 
AOD observations during our case study of 4 August 2007 
are below 1 . 

For our case study, the version 3 total column AOD 
is equal to the version 2 AOD of 0.32. Having studied 
the different potential reasons above for a CALIOP extinc- 
tion underestimation, we speculate that the effects of fixing 


the cloud-contamination bug (reducing the AOD) are being 
counterbalanced by the consequences of a downward exten- 
sion of the layer base altitude for low-lying aerosols (increas- 
ing the AOD) in the version 3 aerosol extinction algorithm. 

In conclusion, this study has helped illustrate potential rea- 
sons for deficiencies in the CALIOP version 2 level 2 aerosol 
extinction product. We hope that our study will improve 
the understanding of the results obtained in previous stud- 
ies that have used CALIOP version 2 data. The next version 
of CALIOP data (version 3) includes corrections to many of 
the factors described above. We emphasize that, as this study 
is not a global validation assessment of the CALIOP version 
2/version 3 aerosol extinction product, it does not provide a 
quantitative guidance regarding the generality and magnitude 
of errors likely to be present in both versions of the CALIOP 
extinction retrieval algorithm. This should be the object of 
future studies. 


Appendix A 

Alternative CALIOP extinction retrieval 


The alternative CALIOP extinction retrieval is based on a 
simple numerical integration technique. Let us define: 

- 0: range of the CALIOP LIDAR (corresponding to alti- 
tude ~705 km) 

- 0— ro: range where there are no aerosols (ro is usually 
equal to r c , calibration range) 

- r : range of the scattering aerosol layer. 

The LIDAR signal, P(r), can be written as follows: 

P(r) = x p(r) x T(0,r) 2 = ^ x p(r) x r(0,r 0 ) 2 x T(r 0 ,r ) 2 (Al) 

where K : the system constant, r: the range, /3(r): the to- 
tal backscatter coefficient profile and T (r) 2 : the atmospheric 
two-way transmittance (i.e. the signal attenuation). 

The total attenuated backscatter coefficient profile, 
can be written as follows: 

p\r)=p(r)xT(r 0 ,r) 2 (A2) 

Substituting Eq. (A2) in Eq. (Al) leads to: 

P(r) = ^ x p'(r ) xT (0,r 0 ) 2 (A3) 

r z 

The combination of Eq. (Al) and Eq. (A3) leads to: 

P(r) xT(0,r) 2 = (i\r) x T (0, r 0 ) 2 (A4) 


with: 


T (0,r) 2 = exp 


-2 J (a a (r')+a m (r'))dr 
0 


(A5) 
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T (0,ro) 2 = exp 


ro 

-2 J a m (r')dr 
0 


(A6) 


P(r) = p a (r) + fi m {r) 


(A7) 


S a (r) = 


Ola 0 ) 

Pair) 


(A8) 


where a a (r) and a m (r) are the aerosol and molecular ex- 
tinction coefficient profiles, /? a 0 r ) and y0 m (?O are the aerosol 
and molecular backscatter coefficient profiles and S a (r) is the 
extinction-to-backscatter lidar ratio. 

Substituting Eqs. (A5-A8) in Eq. (A4) leads to: 


oi a (r) = S a (r) x /f(r)exp 


r 

2 J (a a (r') + a m (r f ))dr 


~P m (r) (A9) 


Let us write Eq. (A9) replacing the integral by a discrete 
sum over different layers i of range r, in the vertical: 


Wa ( r i ) — Sa ( r i ) x 


^'(r/)exp 


" /-I 

2 ^ (a fl (rjfe) + (rjfc)) A r 
_ ife=0 



(A10) 


a a (r/) can then be computed using S a (rQ from the 
collocated HSRL instrument, /3(r/) from the CALIOP 
^532 @1/3 km profile product (cloud- screened with the 
CALIOP cloud @1/3 km product), P m ( r i) and a w (r/) from 
the gridded molecular number density profile data from the 
GEOS-5 analysis product available from the NASA God- 
dard Global Modeling and Assimilation Office (GMAO) and 
knowing A r, the thickness of each layer i on the vertical. 

Here ro is where the aerosol load is assumed to be negligi- 
ble and according to Eq. (A 10), the aerosol extinction coeffi- 
cient ota at range r\ (first layer below ro) can be computed as 
follows: 


Ola {n ) = S a (r 1 ) x (p' (ri ) exp [2a m (r 0 ) Ar ] - p m (ri )) (All) 
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